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ABSTRACT 

We present a spectroscopic analysis of VLT/X-Shooter observations of six O-type stars in 
the low-metallicity (Z ~ 1/7 Z ) galaxies IC 1613, WLM and NGC 3109. The stellar and wind 
parameters of these sources allow us, for the first time, to probe the mass-loss versus metallicity 
dependence of stellar winds below that of the Small Magellanic Cloud (at Z ~ 1/5 Zq) by 
means of a modified wind momentum versus luminosity diagram. The wind strengths that we 
obtain for the objects in WLM and NGC 3109 are unexpectedly high and do not agree with 
theoretical predictions. The objects in IC 1613 tend towards a higher than expected mass- 
loss rate, but remain consistent with predictions within their error bars. We discuss potential 
systematic uncertainties in the mass-loss determinations to explain our results. However, if 
further scrutinization of these findings point towards an intrinsic cause for this unexpected sub- 
SMC mass-loss behavior, implications would include a higher than anticipated number of Wolf- 
Rayet stars and Ib/Ic supernovae in low-metallicity environments, but a reduced number of 
long-duration gamma-ray bursts produced through a single-star evolutionary channel. 
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The evolution of massive stars is greatly af- 
fected by the amount of mass and angular mo- 
mentum lost during their lifetime. Understand- 
ing the mechanisms responsible for these losses is 
a fundamental goal in stellar astrophysics. For 
instance, mass-loss influences the characteristics 
of the supernova explosion with which a mas- 
sive star ends its life, as well as the number of 
potential single-star prog enitors of long-duration 



and driving an outflow ( e.g., Castor et al. 19751 : 



Kudritzki fc Puis 1 120001 IVink et al. Il200ll) . The 



gamma-ray bursts (e.g . Yoon fc Langer 2005; 



Wooslev fc Bloom I 2006). 



An important mass-loss mechanism, at least at 
galactic and Magellanic Cloud metallicities, is the 
transfer of momentum from photons to the at- 
mospheric gas through line interactions, initiating 



1 Based on VLT/X-Shooter observations under program 
085D.0741. 



strength of these radiation-driven winds therefore 
depends on the effective number of absorption 
lines, in particular near the photospheric flux max- 
imum in the ultraviolet. It is dominated by the 
absorption of light through a copious amount of 
metallic ion lines that are present in this wave- 
length region. As a consequence, the mass-loss 
rate (M) of stars hotter than 25000 K is predicted 
to scale with met al licity (Z) as Ma z ° S9±0 - w 
(|Vink et al. Il200l[) . iMokiem et al~l (|2007l) showed 
that this prediction holds for early-type stars in 
the Galaxy (Z = Z Q ), Large Magellanic Cloud 
(LMC; Z = 0.5 Z Q ), and Small Magellanic Cloud 
(SMC; Z — 0.2 Zq), yielding the empirical rela- 
tion M oc Z°- 78±0 - 17 . To date however, no ob- 
servational constraints exist at sub-SMC metallic- 
ities. 
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Table 1 

Properties of the observed stars 



ID 


K.A. 


Dec. 


T 7 
\ 


My 


Spectral 


ix. V . 


lZ+10g(U/riJ 




(J2000) 


(J2000) 






Type 


(km s -1 ) 




IC 1613 
















— A13 


01 05 06.21 


+02 10 44.8 


19.02 


-5.55 


03 V((f)) 


-240 


8.0 


— A15 


01 05 08.74 


+02 10 01.1 


19.35 


-5.11 


09.5 III 


-240 




— Bll 


01 04 43.82 


+02 06 46.1 


18.68 


-5.84 


09.5 I 


-240 




— C9 


01 04 38.63 


+02 09 44.4 


19.02 


-5.44 


08 III((f)) 


-265 




WLM 
















— All 


00 01 59.97 


-15 28 19.2 


18.40 


-6.35 


09.7 la 


-135 




NGC 3109 
















—20 


10 03 03.22 


-26 09 21.4 


19.33 


-6.67 


08 I 


407 


7.8 



Refer ences. — y-band m agnitudes from iBresolin et al. 
WLMl: lEvans et al. I (120071 NGC3109) 



(|2007L IC1613); [Bresolin et al. ] (120061 



For galactic stars more luminous than 10 Lq 
(i.e. above the Humphreys-Davidson limit), 
episodic mass loss is expected to occur dur- 
ing the Luminou s Blue Variable and/or Wolf- 
Rayet phase (e.g., Humphreys fc Davidson Ill994t 



Smith et al. 1 120041 ) . Also, dust/pulsation-driven 
mass loss takes place in the red sup ergiant phase 
of lower luminosity objects (e.g., iHeger et al 



1997t lYoon fc Cantiello I l20ldh . However, if the 



relative contribution of line driving remains siz- 
able at all Z, an important consequence of the 
theory of radiation-driven winds would be that 
low-metallicity massive stars loose less mass and 
angular momentum over their lifetime than their 
high-metallicity counterparts. At sub-SMC mctal- 
licities, single stars with a large initial rotational 
velocity may avoid envelope expansion and the 
associated efficient transfer of angular momen- 
tum from the core to the envelope, thus keeping a 
rapidly spinning core. This makes them potential 
prog enitors of long-duration gamma-ray bursts 
(e.g. IWooslev fc Bloom II2006I ). 

In order to observationally constrain stel- 
lar and wind parameters at sub-SMC metallic- 
ities, we have used the X-Shooter spectrograph 
moun ted on ESO's Very Lar ge Telescope (VLT) 
UT2 (jP'Odorico et al. H2006L Vernet et al. 2011, 
in press) to secure intermediate-resolution spectra 



of some of the most massive stars in three nearby 
dwarf galaxies, IC 1613, WLM and NCG3109, 
each having a metallicity of Z « 1/7 Zq. The 
throughput of the instrument, combined with the 
collecting area of an 8.2m telescope, allows us, for 
the first time, to perform a detailed quantitative 
spectral analysis at a resolution R ~ 6000-9000 
and test the theory of line driving in O-type stars 



at a sub-SMC metallicity (see also iHerrero et al. 
201lh . 

With the aid of pre v ious, l ow- resolution studie s 



(|Bresolin et al. 1 120061 l2007t lEvans et all 120071 ). 



we selected the six visually brightest O-type stars 
in these galaxies (four in IC 1613, one in both 
WLM and NGC 3109). The low line-of-sight ex- 
tinction towards these galaxies and their similar 
metallicity allow us, for this particular study, to 
treat these stars as a group. 

In Section [2] the observations and data reduc- 
tion are described. Section [3] discusses the method 
of analysis and Section |4] presents the results. Fi- 
nally, in Section [5j we discuss implications of our 
findings. 

2. Observations & Data Reduction 

The spectrum of NGC 3109-20 has been ob- 
tained during a X-Shooter GTO run on 2010 April 
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Fig. 1. — Selected regions of the normalized X-Shooter spectra of the observed O-stars. Overplotted in red 
are the best-fit line profiles for the fitted lines. 



23, with an exposure time of 5x900s. In a succes- 
sive GTO run from 2010 September 12 to 16, the 
stars in IC 1613 and WLM were observed, with 
an exposure time of 6x900s for all stars. All ob- 
servations were obtained in nodding mode with 
a nod throw of 5", using a slit width of 0.8" for 
the UVB arm (300-550 nm), 0.9" for the VIS arm 
(550-1020nm), and 0.9" for the NIR arm (1000- 
2500 nm), yielding a spectral resolving power of 
i?=6200, 8800 and 5600, respectively. Conditions 
were clear or photometric with an average seeing 
below 0.8" in the R band. The moon illumination 
fraction was below 0.2. Fundamental properties of 
the observed stars are listed in Table [TJ 

2.1. Data Reduction 

The data have been reduced using the X- 
Shooter pipeline vl. 2 .2 in physical mode l mod e 



( Gold oni et al. ] 120061; iModigliani et al. I l2010h 



The nodding mode usually allows for easy nebular 



subtraction. However, all sources except WLM- 
All show signs of variable strength of the nebular 
emission lines over the length of the nod throw, 
preventing the use of this method. For these stars 
the sky subtraction was done by using suitable re- 
gions of sky close to the object, resulting in a good 
nebular correction for all stars except IC 1613-A15 
(see below). The signal-to-noise ratio (SNR) per 
resolution element ranges from 32 to 67 in the 
UVB arm and from 21 to 29 in the VIS arm. The 
detected flux in the NIR part of the spectrum was 
low (SNR < 5), so we could not use this part of 
the spectrum in our analysis. 

The nebular emission in IC 1613-A15 varies in 
radial velocity along the slit, causing residuals in 
the final spectrum. The regions affected by these 
residuals were removed from the final spectrum 
and not used in the atmosphere fit. 

The resulting ID spectra were normalized to 
the continuum. Figure []] presents the blue spec- 
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trum and Ha region of all six stars. Figure [2] 
shows the profiles of the mass-loss sensitive He II 
4686 and Ha lines in more detail. The spectral 
type of the stars was determined by comparing 
with a spectral atlas of O-type stars (Sana et al. 
2011, in preparation), and for all but one agree 
with the type listed in the literature. The only 
change is that the spectr al type of IC1613-A1 3 is 



He II 4686 



Ha 



20071) to 



refined from 03-4 V((f)) (|Bresolin et al. 
03 V((f)). 



2.2. Oxygen Abundance Determination 

In order to also provide an independent oxygen 
abundance measurement, a relative flux calibra- 
tion was performed for the frames where nebu- 
lar oxygen emission is present (IC 1613-A13 and 
NGC 3109-20), using spectra of photometric stan- 
dard stars taken during each night of observa- 
tions. To derive the oxygen abundance, we use the 
strong line method ( Pagel et al. 1973) with th e 
updated calibration bv lPilvugin fc Thuan I (|2005l ). 
using the [O II] AA 3727,3729, [O III] AA 4959,5007 
and H/3 nebular emission lines. The resulting 
abundances are given in Table [TJ and although 
abundances derived with the strong line method 
can have unce r tainti es up to 0.8 dex (see e.g. 



Bresolin et al. 20091 ) . our measures agree well 
with the literature values of 12 + log(0/H) = 



7.90±0.08 ( IC 1613:lBresolin et al. I liooj. 7.80± 
0.07 fWLM: IUrbaneia et al. 1120081) and 7.76±0.07 
(NGC 3109: lEvans et al. 1 120071) derived from B- 
and A-type supergiants. 

3. Modeling 

The stellar properties and wind characteristics 
of the stars have been determined using an auto- 



mated fitting method developed by lMokiem et al. 
(|2005h . This method combines the non-LTE 



stella r atmosphere model FASTWIND (jPuls et al 



2005 ) with the gene ric fitting algorithm PIKAIA 
( Charbonneau 1995). It allows for a fast and ho- 
mogeneous analysis of our sample of stars using a 
selection of hydrogen, He I, and He II lines. 

The absolute visual magnitude (My) of the 
stars is needed as an input parameter. They 
have been determined using the V magnitudes 
from Table [1] an d a distance and redd ening 
of d 
E(B 



721 kpc (lPietrzvnski et al. 
V) = 0.025 ijSchlegel et al. 



2006) and 
199$) for 




IC 1613-B1 1 



E i.o 



IC 1613-A13 






NGC 3109-20 



WLM-A1 1 




-600-400-200 200 400 600 -600-400-200 200 400 600 

Velocity (km/s) 

Fig. 2. — Line profiles of the mass-loss sensitive 
diagnostic lines He II 4686 and Ha. The best- 
fit line profiles are overplotted with solid lines. 
The dashed profiles indicate models with mass-loss 
rates predicted by radiation-driven wind theory. 



IC 1613, d = 995 kpc and E(B - V) = 0.08 
for WLM (lUrbaneia et al. I l2008l). and d = 1300 
kpc (ISoszyhski et al. 11200(3) and E{B - V) = 0.14 
(jDavidge Ill993l) for NGC 3109. The resulting val- 
ues of My ag ree well with the va lues for their 
spectral type ( Martins et al. 2005). The radial 



velocity (R.V.) has been measured by fitting a 
Gaussian profile to the hydrogen lines, and the 
obtained values are in agreement with the radial 
velocities of the host galaxies. 

The fitting algorithm covers a large parameter 
space, fitting line pr ofiles of the same 11 s pectral 
lines as described by iMokiem et al. ( 2005 ) to de- 
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termine the effective temperature (T e gf), gravity 
(g), mass-ioss rate (M), surface helium abundance 
(Yho), depth-independent microturbulent velocity 
(«tur)) and projected rotational velocity (w ro t sini). 
The bolometric luminosity L is derived by ap- 
plying the bolometric correction to the absolute 
visual magnitude used as input. This luminos- 
ity, together with the obtained temperature, is 
then used to determine the radius (R) of the star. 
The terminal wind velocity (woo) can not be con- 
strained from the optical, but is related to the 
surface es cape velocity ^esc : ypo — ^esc (&t 



Lamers et al. 1995; Kudritzki fc Puis 



2000). The coefficient of the wind velocity struc- 
ture P has been fixed to 0.8 for the dwarfs, 0.9 
for the giants, and 0.95 for the supergiants, con- 
form theoretical predictions (Muijres et al. 2011, 
in preparation). 

Table [2] presents the best- fit values and derived 
properties for the six stars. The values for T c ff 
are higher than those of their G alactic counter- 
parts of the same spectral type ( Martins et al. I 
2005 ), in agreement with their low metallicity 
(|Mokiem et al~ll2004h . IC 1613-A13 shows an un- 
usually high value for Yho • This is likely caused by 
the degeneracy between effective temperature and 
surface helium abundance at high temperatures 
(resulting from the loss of He I as a diagnostic). 
The large uncertainties in the determination of the 
micro-turbulent velocity show that this parameter 
cannot be well constrained with our data. 

The errors on the parameters have been an- 
alyzed by calculating the probability (P = 1 — 
r (x 2 /2, z//2), where F is the incomplete gamma 
function and v the degrees of freedom) for all cal- 
culated models. Because P is very sensitive to the 
value of the % 2 , we normalize all \ 2 values such 
that the best x 2 od is equal to 1, i.e. we assume 
that deviations of the original best x 2 cd from unity 
are induced by under- or overestimated errorbars 
on the normalized flux. This approach is similar 
in spirit to using relative weighting in the \ 2 merit 
function and to propagating the root- mean-square 



of th e fit to scale the error bars fe.g.. iPress et al. 
1986). The uncertainties are obtained by consid- 
ering the range of models which satisfy P > 5%. 
These errors do not take into account uncertain- 
ties in the luminosity. However, as the mass- 
loss rate approximately scales with luminosity as 
M <x L 5 / 4 , typical uncertainties in L do not have 



a significant impact on our conclusions. 

4. Results & Discussion 

The derived mass-loss rates are represented in 
the modified wind-momentum vs. luminosity di- 
agram (WLD, Figure [3]). The modified wind mo- 
mentum is defined as D mom — Mvqo\J ' R/ 'Rq. 
This quantity is ideally suited to study the M{Z) 
relation because it is almost independent of mass, 
and Voo and R are usually relatively well con- 
strained. As the Ha recombination line is essen- 
tially sensitive to the invariant wind-strength pa- 
rameter Q = M I '(i? 3 / 2 ^) that is inferred from 
the spectral analysis, D (for given T e g) scales with 
L, making it less sensitive to uncertainties in the 
luminosity. As Uqo is expected to scale wi th metal- 
licity (woo cx Z 013 ; lLeitherer et al. |[l992l) . and the 
wind-strength parameter Q is invariant, the de- 
rived mass-loss rate is subject to a similar scaling. 
This scaling of M and «oo has been applied to the 
values given in Tabled] 



2001, Z 



Compared to theoretical expectations (jVink et al. 
0.14Z© dashed line in Figure [3j> , all 



stars except IC 1613-A13 tend toward a higher 
than predicted mass-loss rate for their metallic- 
ity. Compared t o the empirical relations found by 
Mokiem et al. I ([20071 ) . our values are reminiscent 
of the values measured for the LMC. For IC 1613- 
A15 we could only obtain an upper limit due to 
the nebular contamination, and IC 1613-B11 and 
C9 could have a low enough mass-loss rate within 
the errors. WLM-A11 and NGC 3109-20 have 
a well defined mass-loss rate which is almost an 
order of magnitude too large for their metallicity. 

4.1. Systemetic Uncertainties in the Mass- 
Loss Determination 

In addition to random uncertainties, several 
sources of systematic uncertainties may affect our 
mass-loss determinations, for instance due to as- 
sumptions we have made. Here we discuss po- 
tential non-intrinsic causes for the high mass-loss 
rates that we derive. 

The adopted values of the flow acceleration pa- 
rameter /? could contribute to the high values of 
D mom . To fit the line profile, the wind density 
vs. velocity profile must be the same, causing a 
higher value of j3 to give a lower mass-loss rate. If 
we underestimated (3 the mass-loss rate will there- 
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Table 2 

Best-fit parameters and derived properties of the observed stars 



ID 


Toff 




logM 




Vtut 


Wrotsin i 


Woo 


logL 


R 




[kK] 


[cm s~ 2 ] 


[Ms y^ 1 ] 




[km s^ 1 ] 


[km s" 1 ] 


[km s^ 1 ] 


[Le] 


[Re] 


IC 1613 




















— A13 


47 fi+ 4 - 73 

Z i'-D_ 4 95 


o 70+O.I3 
°-'°-0.22 


-6 26+ 045 

°- ZD -0.50 


n qn+0.04 
u - oz -0.17 


4+25 


94^ 


1869 


5.78 


11.4 


— A15 


no 7 +3.10 
°°- ' -2.65 


3 76^ 
°-'°-0.39 


-6.36+ 040 


16+ - 14 

U ' 1U -0.10 




34ig 


1971 


5.24 


12.2 


— Bll 


01 O+2.00 
ol -°-2.00 


41+O.24 
J -^ 1 -0.21 


-fi 16+ - 30 

D - 1D -1.30 


u - io -0.07 


12+ 13 


qq+36 
°°-28 


1601 


5.45 


18.1 


— C9 


35-7±i:g 


e.p+0.19 
O.OO_ 24 


-6 26+ 25 

°- zo -0.60 


u - lz -0.04 


15+ 8 


70+28 
' z -28 


1697 


5.43 


13.6 


WLM 




















— All 


29.7l 2 2 . 4 7 5 5 


9 r+0.29 
°- zd -0.19 


_S rc fi +0.20 
°- OD -0.30 


0.11J 11 


8+ 14 




1711 


5.79 


29.8 


NGC 3109 




















—20 


34.2±t™ 


O A O+0.37 


r 41+O.25 
°- 41 -0.35 




16j 


98± 8 7 | 


2049 


5.88 


24.7 



Note. — Arrows indicate upper or lower limits. 



fore be too high. However, even for extreme val- 
ues of (3, the correction will be at most a factor 
of two. The assumption that — 2.6 v OS c intro- 
duces an uncertainty in Z) mom , but as errors in the 
terminal flow velocities are not expe cted to exceed 



20-40% (jGroeneweeen et al. I [1989) . this effect is 



even smaller. 

Another source of uncertainty is the luminosity 
determination, both through extinction and dis- 
tance uncertainties. However, as said before, the 
WLD is not very sensitive to uncertainties in L, 
as an increase in L almost linearly increases the 
value of flmom. This effect is illustrated by the 
arrow in Figure [3J indicating a change of a fac- 
tor of two in distance. Furthermore, since the 
extinction towards all galaxies appears very low, 
and the absolute magnitude of the stars (Table Q} 
agrees well with the v alues for their spectral type 
([Martins et al. Il2005l ). we do not expect a signifi- 
cant error in the luminosity. 

An underestimated metallicity could also be 
part of an explanation for the difference in em- 
pirical and predicted mass-loss rates. However, 
the various studies mentioned in Table Q] all agree 
on an oxygen abundance of 10-15% solar. Our 
own determination of the oxygen abundance for 
IC 1613 and NGC 3109 also agrees with these 
values. Furthermore, earlier studies on the stellar 



popu la tions of our target galax i es (ICole et al. 



1999 : iMinniti fe Ziilstra I Il997t iMinniti et al. 



1999() derived iron abundances of again 10-15% 



solar. 

Finally, multiplicity could play a role. Since 
the observed spectroscopic bi nary fraction in the 
Galaxy is approximat ely 50% ( Mason et al. | [2009; 
Sana fc Eva"ns1l2010h . it is possible that some of 
our sources are binaries. Binarity can influence 
the derived mass-loss rate by dilution or wind- 
wind collisions. In the case of a low-mass com- 
panion, the normalized spectrum will have shal- 
lower lines due to dilution effects. In close bina- 
ries with two massive stars, wind-wind collisions 
can cause emission in He II A4686 and Ha if the 
wind -collision region is isothermal (see, for exam- 
ple, [S^na_eLaL ll200ll) . However, we could not find 
indications for massive companions in our data, 
which would cause higher luminosities, nor of ra- 
dial velocities that differ from those of the parent 
galaxy. 

Because of the large distance to the host galax- 
ies, it is possible that our targets are unresolved 
multiples instead of single stars (e.g., Hartoog et 
al. 2012, in preparation). With a seeing of 0.8", 
the smallest resolved area covers a distance rang- 
ing from 2.8 pc for IC 1613 to 5 pc for NGC 3109, 
and can therefore easily harbor a typical open 
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Fig. 3. — Modified wind-mo mentum vs. lumin osity diagram with our results. The dashed lin e s ind icate 
the theoretical predictions of Vink et al. I ( 200ll ). and the empirical results from Mokiem et al. ( 2007 ) are 
represented by the shaded bars. The arrow illustrates the effect of an uncertainty in the distance on the 
location of points in the diagram. 



cluster. However, the observed luminosities are 
typical of stars of the corresponding spectral type 
( Martins et al. Il2~005() , and we therefore do not see 
any strong indications for multiplicity. 

4.2. Possible Intrinsic Causes for the High 
Mass-Loss Rates 

There are several physical mechanisms that can 
increase the mass-loss rate or make it appear 
higher. A likely cause is the effect of wind clump- 
ing, which we do not take into account in our mod- 
eling. Clumping causes mass-loss indicators based 
on the average squared density (p 2 ), such as Ha, 
to have more emission than would be the case in an 
unclumped wind. This causes the derived mass- 
loss rates to be overestimated by a factor of 1/y/J, 
with / the mean value of t he volume- filling factor 



in the line-forming region (jPuls et al. 1 12008). 



clumping is also unaccounted for (causing the em- 
pirical mass-loss rate to be higher than the the- 
oretical prediction), the rates we find for Z ~ 
0.14 Zq galaxies coincide with LMC values. So 
unless clumping behaves differently at sub-SMC 
metallicities, we would still expect the mass-loss 
rates to be located below the empirical SMC area 
in the WLD. 

The slope and absolute scaling of the WLD is 
also anticipated t o change at very low metallicities 
(jKudritzki 1 120021) . The expected effect, however, 
works in the opposite direction, and the mass- 
loss rates should be even lower. Furthermore, the 
slope-change is predicted to occur at a much lower 
metallicity (Z ks 10" 3 Z q ). 

Fast rotation can also c ause increased and 
asym metric mass-loss (e.g., iMaeder fc Mevnet" 
2000). However, all of the projected rotational 



However, by co mparing our re s ults to the re- 
lations found by iMokiem et al. ( 2007 ). where 



velocities we derive are low (v sini < 100 km s 1 ), 
and given that the probability that by chance all 
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sources have their rotational axis pointed more or 
less in our direction is small, we do not expect this 
effect to play a major role in our sample. 

There could be other processes affecting the 
mass-loss process, possibly related to pulsa- 
tions and magnetic fields, which both have 



been detected 
19991: 



O-type 



Donati &: Landstreet 



st ars ( e.g. 
20Q9h . 



Henrichs 



Pulsations 

are caused by the high sensitivity of opacity to 
temperature in the high-temperature low-density 
regions in O (super)giants, givin g rise to k- 
pulsations (e.g. Ilglesias et al. 1 1 19921 ) . This is in- 
duced by the sudden appearance of a large number 
of same-shell transition iron lines, and is therefore 
dependent on the iron abundance. Thus it is ex- 
pected that pulsations are less im portant at low 



mctallicitics (Bara ffe et al. I 2001) 



The increased opacity described above can give 
rise to small convect ive reg ions in the stellar enve- 
T200E 



Cantiello fc Braithwaite 



lope (ICantiello et al. 
(|2011l ) showed that these convective regions can 
produce magnetic hot spots through the dynamo 
effect, which can be strong enough to influence 
the wind and possibly play a role in the wind 
clumping. This effect is again metallicity depen- 
dent, and is expected to play a lesser role and 
eventually disappear w ith decreasing metallicity 



(Cantiello et al. 



5. Conclusions and implications 

In this Letter, we have pointed out a discrep- 
ancy between the observed and predicted mass- 
loss rates from massive stars in low-metallicity en- 
vironments, and discussed possible explanations. 

A potential violation of the expected metallic- 
ity scaling of the radiation-driven mass-loss rate 
at Z ~ 1/7 Zq, resulting in higher than expected 
mass-loss rates, would have far-reaching implica- 
tions. Firstly, one expects \ow-Z O-type stars to 
suffer more from spin-down through angular mo- 
mentum loss in the stellar wind. Consequently, 
the rotational mixing efficiency may be reduced, 
leading, for instance, to a more mo dest nitrogen 



enric hment than currently thought (jBrott et al. 
20111) . 



Secondly, one anticipates the single O star pop- 
ulation in low metallicity environments to produce 
more observationally identifiable Wolf-Rayet stars, 
being the successors of O-type stars having their 



outer envelope stripped by mass-loss, and there- 
fore increased number of lb and, potentially, Ic su- 
pernovae. The single-star channel would, however, 
produce less progenitors of long-duration gamma- 
ray bursts, as the stars loose more angular mo- 
mentum by their outflow. 

Finally, if the larger than expected wind 
strength at Z ~ 1/7 Zq persists to extremely low 
metallicities (though at present the driving mech- 
anism is unknown), stellar winds will impact the 
evolution of massive Population III stars and the 
chemical enrichment of the intergalactic medium 
of the early Universe. 

We thank the referee, dr. Kudritzki, for his 
useful comments, dr. Martayan for his support 
during the observations, and O. Hartoog for the 
interesting discussions. 

Facilities: VLT:Kueyen (X-Shooter) 
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